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ABSTRACT 

The  Center  Manifold  Theorem  is  applied  to  the  local  feedback  stabilization  of  nonlinear 
systems  in  critical  cases.  The  paper  addresses  two  particular  critical  cases,  for  which  the 
system  linearization  at  the  equilibrium  point  of  interest  is  assumed  to  possess  either  a  simple 
zero  eigenvalue  or  a  complex  conjugate  pair  of  simple,  pure  imaginary  eigenvalues.  In  either 
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form,  i.e.,  before  reduction  to  the  center  manifold.  Moreover,  the  formulation  given  in  this 
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and  those  of  an  alternative  approach. 
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1.  INTRODUCTION 


Recently,  center  manifold  reduction  has  been  employed  in  nonlinear  stabilization,  re¬ 
sulting  in  stabilizing  control  law  designs  for  various  classes  of  nonlinear  systems  in  the 
so-called  “critical  cases.”  Critical  cases  occur  when  the  linearized  system  at  an  equilibrium 
point  has  at  least  one  eigenvalue  on  the  imaginary  axis,  with  the  remaining  eigenvalues  in 
the  open  left  half  of  the  complex  plane. 

Aeyels  [1],  who  initiated  application  of  the  center  manifold  reduction  in  nonlinear  sta¬ 
bilization,  investigated  the  existence  of  smooth  stabilizing  feedback  control  laws  for  a  class 
of  third-order  nonlinear  systems  for  which  the  linearized  model  possesses  an  uncontrollable 
pair  of  pure  imaginary  eigenvalues.  Behtash  and  Sastry  [10]  used  the  same  approach  to 
study  stabilization  for  nonlinear  systems  whose  linearized  model  has  two  distinct  pairs  of 
complex  conjugate  pure  imaginary  eigenvalues,  or  a  double  pole  at  the  origin,  or  a  pole 
at  the  origin  and  a  complex  conjugate  pair  of  pure  imaginary  eigenvalues.  In  [10],  the 
design  was  undertaken  for  the  reduced  system  on  the  center  manifold  using  normal  form 
calculations,  and  for  simplicity,  a  scalar  stable  mode  was  assumed.  Generally,  there  is  a 
need  for  considering  cases  with  any  finite  number  of  stable  modes.  Moreover,  it  is  desirable 
to  express  the  control  laws  directly  in  terms  of  the  original  model  rather  than  in  terms  of 
transformed  versions. 

A  main  goal  of  this  paper  is  to  derive  such  stabilizing  control  algorithms  for  general 
nonlinear  systems  in  critical  cases.  The  development  focuses  on  general  nonlinear  systems  in 
two  specific  critical  cases.  In  the  first  critical  case  of  interest  here,  a  simple  zero  eigenvalue 
occurs,  while  in  the  second  case  a  pair  of  pure  imaginary  eigenvalues  occurs.  In  either  case, 
the  critical  eigenvalues  of  the  linearized  model  need  not  be  controllable.  The  feedback 
laws  obtained  include  purely  linear  state  feedbacks,  purely  nonlinear  state  feedbacks  and 
feedback  control  laws  containing  both  linear  and  nonlinear  terms  in  the  state.  Results  of 
this  paper  are  compared  with  those  of  [6],  [7]. 

2.  PRELIMINARIES 

Consider  a  class  of  nonlinear  autonomous  systems 

V  =  Anri  +  Al2£  +  F(r},()  (la) 
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f  =  Anr}  +  A22C  +  G(i],  £),  (lb) 

where  i]  £  IR11 ,  £  €  IRm .  In  (1),  A,y  for  i,j  —  1,2  are  constant  matrices,  and  the  functions 
F,  G  are  sufficiently  smooth,  with  their  values  and  first  derivatives  vanishing  at  the  origin. 
If  A12  and  A21  vanish,  the  matrix  An  has  all  its  eigenvalues  on  the  imaginary  axis,  and  A22 
is  Hurwitz,  then  the  Center  Manifold  Theorem  asserts  the  existence  of  a  locally  invariant 
manifold  for  (1)  near  the  origin.  This  manifold  is  given  by  the  graph  of  a  function  £  =  h(rj). 

In  applying  the  Center  Manifold  Theorem  to  feedback  stabilization  problems,  it  is 
convenient  to  give  a  restatement  of  the  theorem  in  a  way  that  does  not  require  vanishing 
of  the  “linear  coupling”  matrices  A12  and  A21.  This  is  especially  true  when  the  feedback 
is  allowed  to  possess  linear  terms.  For  the  purposes  of  this  paper,  a  restatement  allowing 
nonzero  A21  but  with  A12  =  0  suffices.  A  linear  transformation  of  variables  is  now  employed 
to  achieve  this.  Consider  the  equation 

AM  +  MB  =  C,  (2) 

where  A  £  (PmXm,  B  £  <SnXn  and  M,  C  £  <£mXn .  For  n  =  m  and  B  =  Ar,  Eq.  (2)  is  a 
Liapunov  matrix  equation  [5].  Let  T  denote  the  linear  operator 

T  :  M  >->  AM  +  MB  (3) 

for  M  £<SmXn. 

The  following  result  is  a  direct  generalization  of  [5,  Theorem  F-l  and  Corollary  F-lcr.]. 

Theorem  1.  Let  n,m  be  positive  integers.  If  the  sum  of  any  eigenvalue  of  A  and  any 
eigenvalue  of  B  is  nonzero,  then  the  linear  matrix  equation  (2)  has  a  unique  solution  for 
matrix  M. 

We  now  apply  the  Center  Manifold  Theorem  to  the  stability  analysis  of  (1)  for  the 
case  in  which  A12  =  0,  with  A21  not  necessarily  zero.  Let  A22  be  Hurwitz  and  Au  have  all 
its  eigenvalues  on  the  imaginary  axis.  By  Theorem  1,  the  equation 

A22E  —  EAu  A  A21  =0  (4) 

has  a  unique  solution  for  the  m  x  n  matrix  E.  Letting  v  :=  £  —  Erj ,  we  can  rewrite  system 
(1)  as 

V  =  An  77  +  F(r),v  +  Erj)  (  5  a) 

v  =  A2 2V  +  G(rj,  v  -f  Eif)  —  E  ■  F(rj ,  v  +  Erj).  (5b) 
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The  Center  Manifold  Theorem  for  (1)  can  now  be  restated  as  follows: 

Lemma  1.  Assume  A\2  =  0,  A22  is  Hurwitz,  and  all  eigenvalues  of  An  have  zero  real  parts. 
Then  the  origin  of  (1)  is  asymptotically  stable  (unstable)  if  the  origin  is  asymptotically 
stable  (unstable)  for  the  reduced  model 

V  =  Anrj  +  F(r),  h(rj)  +  Er]),  (6) 

where  h  satisfies  the  partial  differential  equation 
Dh(rj){Anr]  +  F(i],  h(rj )  +  Er])} 

=  A22h(r])  +  G(rj,  h(r})  +  Erj )  -  E  ■  F(r 7,  h(rj)  +  Er/)  (7) 

with  E  the  solution  of  Eq.  (4)  and  boundary  conditions:  h( 0)  =  0  and  Dh( 0)  =  0. 

We  employ  Taylor  series  expansions  in  the  development  below,  using  multilinear  func¬ 
tion  notation  for  the  terms  in  these  expansions.  The  definition  of  multilinear  function  is 
recalled  as  follows. 

Definition  1.  (e.g.,  [9])  Let  V\,  V2, . . . ,  14  and  W  be  vector  spaces  over  the  same  field.  A 
map  xp  :  V\  X  V2  X  . . .  x  Vk  — »  W  is  multilinear  (or  ^-linear)  if  it  is  linear  in  each  of  its 
arguments.  That  is,  for  any  rq,  v,  €  Vi,  i  =  1, . . . ,  &,  and  for  any  scalars  a,  a,  we  have 

-0(ui  ,...,avi  +avi,...,v  k)  =  atp(vi , . . . ,  vif . . . ,  vk) 

+  aif>(vi,...,vi,...,vk).  (8) 

The  integer  k  is  the  degree  of  the  multilinear  function  ip. 

The  next  definition  deals  with  the  special  case  in  which  Vf  =  V2  =  •  •  •  =  Vk  -=  V. 

Definition  2.  [9]  A  A:-linear  function  ip  :  V  x  V  x  ...  x  V  —>  W  is  symmetric  if  the  vector 
ipiy i,V2, . . . ,  vk)  is  invariant  under  arbitrary  permutations  of  the  argument  vectors  u,-.  A 
function  <p  :  IRn  — >  IB"1  is  homogeneous  of  degree  k  ( k  an  integer),  if  for  each  scalar  a, 
(pipirf)  =  Oik<p(r])  for  all  rj  6  IRn . 

Note  that,  in  the  sequel  prime  denotes  the  transpose  of  both  vector  and  matrix  and  I 
denotes  the  identity  matrix. 
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3.  GENERAL  FRAMEWORK 


Consider  a  nonlinear  control  system 

V  =  Aur]  +  biu  +  F(r},£),  (9a) 

i  =  ^4.22^  +  b2u  +  G(rj,  £),  (96) 

where  ?7,£  are  real  vectors,  and  a  preliminary  block  diagonalization  has  been  applied  to 
remove  any  linear  coupling  term  in  the  dynamics  between  r]  and  £.  For  simplicity,  u  is 
supposed  to  be  a  scalar  control.  It  is  not  difficult  to  extend  the  study  to  the  case  in  which 

the  input  is  a  vector  control.  In  the  following,  we  apply  the  center  manifold  result  in 

Lemma  1  to  design  stabilizing  control  laws  for  (9)  for  which  all  eigenvalues  of  An  lie  on 
the  imaginary  axis. 

Let  us  first  consider  the  case  in  which  b\  is  nonzero.  In  the  simple  critical  cases,  where 
An  is  the  scalar  0  or  is  a  2  x  2  matrix  with  a  pair  of  pure  imaginary  eigenvalues,  linear 
theory  will  imply  the  existence  of  a  linear  stabilizing  feedback  control  for  (9).  Consider 
next  the  existence  of  a  purely  nonlinear  smooth  feedback  (i.e. ,  one  with  vanishing  linear 
part). 

Since  now  we  focus  on  purely  nonlinear  stabilizing  controllers,  system  (9)  retains  the 
linear  decoupling  property  upon  control.  Thus,  if  A22  is  stable,  then  according  to  center 
manifold  theorem  (e.g.,  [3],  [8])  there  is  a  locally  invariant  manifold  £  =  /i(?/)  for  (9). 
Furthermore,  h  satisfies 

Dh(rj){Ani)  +  biu(rj,  h(rj))  +  F(r /,  h(rj))} 

=A22 h(r])  +  b2u(rj,  h(rj))  +  G(i],  h(r]))  (10) 

with  boundary  conditions  h(0)  =  0  and  Dh(0)  =  0.  Then,  we  seek  a  purely  nonlinear 
stabilizing  feedback  control  law  by  using  stability  conditions  for  the  reduced  model 

r)  =  AuT]  +  biu(r],  h(rj))  +  F(rj,  h(r])).  (11) 

Note  that,  for  the  case  in  which  A22  is  not  stable,  a  linear  state  feedback  Ar>£  is  needed  to 
first  stabilize  A22  +  b2K2- 
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Next,  consider  the  case  of  b\  =  0  and  assume  the  feedback  control  to  be  of  the  form 


u(r1,0  =  K1n  +  K1(  +  U(rl,0,  (12) 

where  £/(-,•)  is  a  smooth,  purely  nonlinear  function  whose  first  derivatives  vanish  at  the 
origin.  Rewrite  the  system  dynamics  (9)  as 

V  =  Anrj  +  F(ri,{),  (13) 

£  =  b2hii]  +  (d-22  +  62 K2)£  +  b2U(rf,  £)  +  G(?/,  £).  (14) 

From  Eq.  (14),  the  feedback  has  given  rise  to  a  linear  coupling  term  between  77  and 
£  in  the  dynamics.  As  discussed  in  preceding  section,  there  is  a  constant  matrix  E  such 
that,  with  v  :=  £  —  Ejj,  the  transformed  version  of  the  control  system  (13)-(14)  is  in  block 
diagonal  form.  Here,  E  is  the  (unique)  solution  of  the  Liapunov-like  equation 

b2E\  +  (A22  4~  b2K2)E  —  EA\  1  =  0.  (15) 

We  assume  that  A22  +  b2K2  is  stable.  Moreover,  since  all  the  eigenvalues  of  An  lie  on 
the  imaginary  axis,  then  Theorem  1  guarantees  existence  of  a  solution  E  to  Eq.  (15).  The 
transformed  dynamics  in  the  states  77  and  £  is  then 

r)  =  Anri  +  F(ri,u  +  Er/),  (16ci) 

v  —  (A22  +  b2K2)v  +  b2U(r ],  v  +  Er /)  +  G(i 7,  u  +  Erj).  (16 b) 

Eq.  (16)  has  a  center  manifold  given  by  the  graph  of  a  function  v  =  h(rj),  where  h 
satisfies 

Dh(rj){Anr)  +  F(i 7,  h(rj)  +  Er])}  =  (A22  +  b2K2)h(r] ) 

+  b2U(r],  h(rj)  +  Er/)  +  <2(77,  h(r])  +  Ei])  (17) 

with  boundary  conditions  h( 0)  =  0  and  Dh( 0)  =  0. 

Lemma  1  implies  asymptotic  stability  of  the  origin  for  (16)  if  the  control  gains  Kx ,  K2 
and  the  nonlinear  function  U  are  chosen  such  that  (i)  A22  +  b2K2  is  Hurwitz,  and  (ii)  the 
origin  of  reduced  model  (16a)  with  u  =  h(r) )  is  asymptotically  stable. 
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We  now  proceed  to  consider  two  special  cases  in  which  the  system  has  only  simple 
critical  modes  (i.e.,  one  zero  eigenvalue  or  a  pair  of  pure  imaginary  eigenvalues)  and  the 
rest  of  the  eigenvalues  are  stabilizable. 

4,  ONE  ZERO  EIGENVALUE 

In  this  section,  we  first  consider  stability  conditions  for  scalar  systems  with  a  zero 
eigenvalue.  These  conditions  are  then  employed  in  the  design  of  stabilizing  control  laws  for 
higher  order  systems  with  a  simple  zero  eigenvalue. 

Consider  a  scalar  real  nonlinear  system 

x  =  dx2  +  ex 3  +  •  •  • .  (IS) 

Stability  conditions  for  system  (IS)  are  given  next. 

Lemma  2.  The  origin  is  asymptotically  stable  for  system  (18)  if  d  =  0  and  e  <  0.  The 
origin  is  unstable  for  (18)  if  d  ^  0. 

Now  consider  Eq.  (9),  with  the  scalar  x  replacing  the  critical  state  7],  and  with 
f(x,0  :=F(x,0 

=  fxxX  -f-  xfx£^  T  £  ./*£££  +  fxxxX  T  X  fxxt^, 

+  xt'fzKZ  +  fedt,te)  +  0(  ||(*, Oil4),  (19) 

G(x,  0  =x2Gxx  +  xGx +  G^(£,  £)  +  x*GXXx 
+  x2Gxx +  xGXtf((, £)  +  G^( £,  GO 

+  o(||(x,0!!4).  (20) 

The  coefficients  in  the  Taylor  series  expansions  (19)-(20)  are  either  constants  or  symmetric 
multilinear  functions  of  their  arguments.  For  instance,  f^{000  and  G^(00  denote  a 
symmetric  trilinear  scalar  function  and  a  bilinear  vector  function  of  £,  respectively. 

In  the  remainder  of  this  section,  stabilizing  control  laws  will  be  obtained  for  system 
(9)  under  one  or  the  other  of  the  following  two  hypotheses. 

Hypothesis  1A.  The  matrix  An  =  0  is  a  scalar  and  b\  ^  0. 

Hypothesis  IB.  The  matrix  An  =  0  is  a  scalar  and  b\  =  0. 
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4.1.  The  case  h  0 


In  this  subsection,  we  consider  the  case  in  which  Hypothesis  1A  holds.  The  control  law 
is  taken  to  be  purely  nonlinear.  Existence  of  a  linear  stabilizing  feedback  for  this  case  is 
evident.  Nonlinear  feedback  controllers  are  none  the  less  desirable  in  certain  applications. 
We  assume  A22  is  stable  and  the  scalar  control  input  is  of  the  form 

u(x,0  =U(x,Q 

. — UxxX  “t"  ^^x^f  T  £  T  ^111^ 

+  x2uxx£  +  x^'ux^(  +  0-  (21) 

According  to  center  manifold  theorem,  the  stability  of  the  origin  for  (9)  coincides  with 
the  stability  of  the  origin  for  the  reduced  model 

x  —  biu(x,h(x ))  +  F(x,h.(x)).  (22) 

Here,  h  solves  Eq.  (10)  with  rj  replaced  by  x  and  with  boundary  conditions  h( 0)  =  0  and 
Dh( 0)  =  0.  Indeed,  solving  (10)  we  have 

h(x)  =  x2hxx  +  0(|:r|3),  (23) 

where 

hXx  =  ^22  2uxx  T  Gxx).  (24) 

From  Lemma  2,  we  now  have 

Lemma  3.  Let  A22  be  stable.  Under  Hypothesis  1A,  the  origin  is  asymptotically  stable 
for  (9)  if  fxx  T  b\uxx  =  0  and  fxxx  T  b\uxxx  —  (fx£  T  b\uX(f)A22  (Om  T  b2'uxxs)  <C  0. 

It  is  obvious  from  Lemma  3  that  a  purely  quadratic  stabilizing  control  law  exists. 

Corollary  1.  Assume  that  A22  is  stable.  Under  Hypothesis  1A,  the  origin  of  (9)  is 
asymptotically  stabilizable  by  a  purely  quadratic  feedback  of  the  form  u  =  uxxx 2  -f  xux^ 

if  A22  Gxx  ^  0. 

Furthermore,  below  we  have  a  purely  cubic  stabilizing  controller  for  system  (9)  when 

fxx  =  0. 
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Corollary  2.  Assume  that  A22  is  stable  and  fxx  =  0.  Under  Hypothesis  1A,  the  origin  of 
(9)  is  asymptotically  stabilizable  by  a  purely  cubic  feedback  of  the  form  u  =  uxxx. r3. 

For  the  case  in  which  A 2 2  is  not  stable,  a  linear  feedback  K2^  is  needed  to  guarantee 
the  existence  of  a  locally  invariant  manifold.  Then  the  design  of  stabilizing  control  laws 
proposed  in  Lemma  3  and  Corollaries  1  and  2  can  be  applied  directly. 

4.2.  The  case  61  =  0 

Next,  we  consider  the  case  in  which  Hypothesis  IB  holds  and  consider  feedback  control 
has  the  form  as 

u(x,Q  =  k1x  +  K2Z  +  U(x,Z)  (25) 

with  k\  a  scalar  control  gain  and  the  nonlinear  control  function  U  as  in  (21). 

From  Section  3,  the  stability  of  control  system  (9)  in  this  critical  case  coincides  with 
the  stability  of  the  reduced  model 

x  -  /(.x,  h(x)  +  Ex),  (26) 

where  E  and  h(-)  solve  Eqs.  (15)  and  (17),  respectively,  under  the  conditions:  (A22  +  b2 K2 ) 
is  stable  and  r]  and  K\  are  substituted  by  x  and  k\ ,  respectively. 

As  above,  h  is  as  in  (23).  We  assume  that  A22  +  b2I(2  is  stable.  Solving  Eqs.  (15)  and 
(17),  we  have 

E  =  — (A22  +  b2I\2)~1b2k1,  and  (27) 

hxx  =  (A22  +  b2K2)  1  {[fxx  +  fx$E  +  E'  f^E]E  —  [b2uxx 
AGXX  +  (b2ux£  +  GX{)E  +  b2E'u^E  +  G^(E,  E1)]}.  (28) 

The  reduced  model  (26)  is  then  given  by 

£  —  {fxx  +  fx^E  +  E'  f^E}x2  +  {fx^hxx  +  2  E'f^hxx  +  fxxx 
+  fxx^E  +  E'fxUE  +  f^(E,  E,  E)}x3  +  0(|x|4).  (29) 

Note  that  E  and  h  as  given  in  (27)  and  (28)  depend  on  the  control  u.  Using  Lemma 
2,  we  have  the  following  stabilization  result  for  control  system  (9). 


9 


Lemma  4.  Let  the  control  input  u  be  of  the  form  as  in  (25).  Then  under  Hypothesis  IB, 
the  origin  of  the  closed- loop  system  (9)  is  asymptotically  stable  if  A22  +  b2 K‘z  is  stable  and 
following  two  conditions  hold: 


fxx  +  fx$E  +  E'  f^E  —  0,  and 

(30) 

fx^hxx  -f-  2 E  f^hxx  fxxx  +  fxx(iE 

+  E'  fx^E  +  f^(E,  E,  E)  <  0. 

(31) 

where  E  is  as  in  (27)  and 

hxx  =  —  (-4-22  +  ^2-ft  2  )  1{i>2 Uxx  +  Gxx 

+  (b2ux^  +  Gxz)E  +  b2E'u^E  +  Gtf(E,  E)}.  (32) 


From  Eqs.  (27)  and  (32),  and  the  fact  that  A22  is  invertible,  we  have  E  =  0  and  hxx  — 
—A22GXX  for  the  uncontrolled  system.  The  next  stability  criterion  for  the  uncontrolled 
version  of  system  (9)  follows  readily  from  Lemma  4. 

Corollary  3.  Suppose  Hypothesis  IB  holds.  Then  the  origin  is  asymptotically  stable  for 
(9)  (with  u  =  0)  if  A22  is  stable,  fxx  =  0  and  fxxx  —  fx^A22Gxx  <  0. 

In  the  rest  of  this  subsection,  we  assume  that  the  stability  conditions  given  in  Corollary 
3  do  not  hold,  and  seek  stabilizing  control  laws  for  system  (9). 

Linear  stabilizing  control  laws  follow  readily  from  Lemma  4,  and  are  as  given  next. 

Proposition  1.  Suppose  Hypothesis  IB  holds  and  let  M  :=  (A2 2  +  b2I(2)~1.  Then  there 
is  a  purely  linear  feedback  which  asymptotically  stabilizes  the  origin  of  (9)  if  there  exist 
feedback  gains  k\  and  K2  for  which  (A22  +  b2K2 )  is  stable, 

fxx  -  h  fx^Mb2  +  k\b'2M' f#Mb2  =  0,  and  (33) 

fxxx  —  fx^MGxx  +  k\  { fx^MGx £  +  2 G'XXM'  fez  —  fxx{\Mb2 

+  kl{b'2 M'fxUMb2  -  fx^M Gtf(Mb2 ,  M&2 ) 

-  2b'2M' f^MGx$Mb2 }  -  k31{f^(Mb2,Mb2,Mb2) 

-  2b2M' f^MG^(Mb2,Mb2)}  <  0.  (34) 
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Remark  1.  The  linear  stabilizing  control  rule  proposed  in  Proposition  1  is  a  composite- 
type  controller  design.  First,  the  feedback  gain  AT2  is  chosen  to  stabilize  state  £.  Then  the 
remaining  feedback  gain  k\  is  selected  to  satisfy  the  conditions  (32)  and  (33)  based  on  the 
chosen  gain  AT2  • 

Since  k\  is  a  scalar,  conditions  (32)  and  (33)  do  not  necessarily  hold  for  any  given  Ah . 
Thus,  a  stabilizing  linear  feedback  does  not  always  follow  from  Corollary  2.  A  special  result, 
for  the  case  in  which  the  non-critical  state  £  is  a  scalar,  is  given  below  to  demonstrate  such 
a  design  is  not  vacuous.  Note  that  £?{$(£,£)  :=  G^£2  in  the  next  corollary. 

Corollary  4.  Suppose  the  non-critical  state  £  is  a  scalar  and  Hypothesis  IB  holds.  Then 
there  is  a  purely  linear  feedback  which  asymptotically  stabilizes  the  origin  of  (9)  if  either 
of  the  following  conditions  holds: 

(i)  ftt  =  OJx*  +  0  and  fxiGxx  -  fxxGxi  +  ±-G^f2xx  <  0. 

(ii)  ±  0,  -  4 fxgfe  >  0  and  either  Gxx  +  GxiE+  +  G^(E+)2  <  0  or  Gxx  +  G^E~  + 

G^(E~)2  >  0,  where 


E  —  g fxt  ^  yj fx£  ~  ^fxxfzs}-  (35) 

According  to  the  stability  conditions  given  in  Lemma  4,  the  cubic  terms  of  both  the 
function  G  and  the  control  input  u  do  not  contribute  to  the  stability  criteria  of  system  (9). 
A  general  linear-plus- quadratic  feedback  control  law  can  then  be  abstracted  as 

u(x ,  £)  =  k\x  +  Af2£  +  uxxx2  +  xux e£  +  (36) 

while  the  control  gains  satisfying  the  conditions  of  Lemma  4. 

As  implied  by  Lemmas  1  and  2  and  the  discussions  above,  we  have  the  next  result. 

Lemma  5.  Suppose  A22  is  stable  and  Hypothesis  IB  holds.  Then  there  exists  no  purely 
quadratic  feedback  stabilizer  for  the  origin  of  system  (9)  if  fxx  ^  0.  However,  the  origin  of 
(9)  is  asymptotically  stabilizable  by  a  purely  quadratic  feedback  of  the  form  u  —  uxxx 2  if 
fxx  =  0  and  fx^A^b2  ^  0. 
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Note  that  the  stabilization  results  given  in  Corollaries  1  and  2  and  Lemma  5  agree 
with  those  obtained  in  [7]. 

5.  PAIR  OF  PURE  IMAGINARY  EIGENVALUES 

In  this  section,  we  consider  system  (9)  in  which  An  has  a  pair  of  pure  imaginary 
eigenvalues.  Specifically,  we  take  An  to  be  of  the  form  (38)  below. 

First,  however,  consider  the  stability  of  a  planar  system 

V  =  An  V  +  Q(v,n)  +  C(rj,rj,rj)  +  •  ■  ■ ,  (37) 


where  ?/  =  (x,y)' ,  and 


An  = 


(38) 


with  17x^2  >  0-  Without  loss  of  generality,  we  may  express  the  quadratic  and  cubic  terms 
in  Eq.  (37)  in  the  form 


Q(v,q) 

C(v,V>v) 


qu  x2  +  q12xy  +  quy2  \ 
721  x2  +  q22xy  +  723 V2  ) 

? 

(39) 

cnx 3  +  cX2x2y  +  c13xy2 
c2xx3  +  c22x2y  +  c2Zxy2 

+  Cl4!/3  \ 

+  C2473  )  ’ 

(40) 

respectively.  Note  the  linearization  of  (37)  at  the  origin  has  the  pair  of  pure  imaginary 
eigenvalues  Ai\fAAh,  where  i  —  \/— T. 

Applying  a  general  stability  criterion  for  planar  systems  undergoing  Hopf  bifurcation 
(see,  e.g.,  [8]),  we  find  that  a  sufficient  condition  for  the  stability  of  the  origin  for  (37)  is: 


-  912(^91.  +  §1«13)  + 

2^2  ,  ^2  1  .  , 
- q2~713723  +  3(cjl  +  3jY"Cl3  +  3C22  +  0~C24^  < 


(41) 


In  the  following,  we  apply  the  stability  criterion  (41)  to  the  design  of  stabilizing  control 
laws  for  the  more  general  (nonplanar)  system  (9)  in  which  both  rj  =  (x,  y)'  and  /q  :  = 
(hi, £>12)'  are  two-dimensional  vectors,  and  F(y,()  =  (f(x,  y,  £),  g(x,  y,  £))' . 
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Results  obtained  in  this  section  will  apply  under  one  or  the  other  of  the  following  two 
hypotheses. 

Hypothesis  2A.  The  matrix  An  (appearing  in  (9))  is  a  2  x  2  matrix  of  the  form  (38) 
above,  and  b\  ^  0. 

Hypothesis  2B.  The  matrix  An  (appearing  in  (9))  is  a  2  x  2  matrix  of  the  form  (38) 
above,  and  b\  =  0. 

5.1.  The  case  6i  ^  0 

First,  we  consider  the  case  in  which  at  least  one  of  bn  and  b\ 2  is  nonzero.  Although  this 
assumption  guarantees  the  controllability  of  the  subsystem  (9),  here  we  consider  only  purely 
nonlinear  control  laws.  Assume  that  A22  is  stable  and  the  control  input  u  =  U(x,y ,£)  is 
a  smooth,  purely  nonlinear  function.  From  Section  3,  the  stability  of  the  origin  of  (9)  now 
coincides  with  the  stability  of  the  origin  of  the  reduced  model 

x  =  +  bnU(x,y,h(x,y))  +  f{x,  y,  h(x,  y))  (42) 

y  =  -&2X  +  b12U(x,y,h(x,y))  +  g(x,y,  h(x,y)).  (43) 

Here,  h  solves  Eq.  (10)  with  y  replaced  by  {x,y)'  and  with  boundary  conditions  h( 0)  =  0 
and  Dh( 0)  =  0.  Indeed,  h  takes  the  form 

h{x,y)  =x2hxx  +  xyhxy  +  y2hyy  +  0(||(®,y)||3),  (44) 

where  hxx,hxy,hyy  are  constant  vectors. 

In  the  following,  we  restrict  the  nonlinear  control  function  U  to  be  a  function  of  x  and 
y  only,  as  follows: 

y^  0  — ^xxX  ~b  TJ>xyXy  T  ^yyV  T  UxxxX 

T  uxxyx  y  T  uxyyXy  T  uyyyV  ■  (45) 

A  stability  criterion  for  the  control  system  (9)  in  this  case  is  given  next. 

Lemma  6.  Suppose  A2 2  is  stable  and  Hypothesis  2A  holds.  Then  the  origin  is  asymptot¬ 
ically  stable  for  (9)  if 

(b\2Uxy  T  9xy){  ( ^12^xx  T  9xx)  T  7T”  (  2  ^yy  T  9yy)} 

o  62 
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where 


(ybl2'UXy  -f-  fxy){  Q  (bllUxx  4~  /n)  "1"  ^2  (.bllUyy  "l-  fyy)} 

2  2^2 
4"  ^  (6llUxx  4“  fxx){bl 2uxx  4“  9xx)  q2  (^11  uyy  4"  fyy){.b\2Uyy  4“  9yy) 

4“  3{&l 1 UXXx  4-  fxxx  4“  fx^hxx  4”  of-v  ( b\\UXyy  4“  fxyy  4“  fx^hyy  4“  fy^hxy) 

Oi£l 

4-  b12uXXy  +  g  xxy  4-  gx£hxy  4-  gi i^h  XX  ) 

0 


+  Q^12  uyyy  +  9yyy  +  9y(b-yy)}  <  0, 


(46) 


hXy  —{^-22  T  1  {2Q,2{Uyyb2  4"  Gyy) 

2£li(uxxb2  4“  G u)  -4-22  (^£1/^2  4"  Gb^)}, 

hxx  =  ■‘^■22  (^££^2  4“  GXx  4“  G,2hxy')i 

b'yy  —  -^22  ( y y b‘2.  4“  Gyy  Q\hXy^) . 


(47) 

(48) 

(49) 


It  is  observed  from  Lemma  6,  generically  there  exists  a  quadratic-plus-cubic  feedback 
stabilizer  for  system  (9).  In  addition,  a  purely  quadratic  state  feedback  stabilizing  control 
law  and  a  purely  cubic  state  feedback  stabilizing  control  law  follow  readily  from  Lemma  C 
as  given  in  the  next  two  corollaries. 


Corollary  5.  Let  A22  be  stable  and  Hypothesis  2A  hold.  Then  the  origin  of  system  (9)  is 
stabilizable  by  a  purely  quadratic  state  feedback  of  the  form  u  —  uxyxy  if 


^>12^Q,1^yy  fxx)~^~  {\2^xx  ^2  fyy } 

—  g  {G2(2gyt  —  8 fxt)A22  4-  +  9x{\(A22  +  4fbfl2-0  1  A22b2  ^  0- 


'50) 


Corollary  6.  Let  A22  be  stable  and  Hypothesis  2A  hold.  Then  the  origin  of  system  (9)  is 
stabilizable  by  a  purely  cubic  state  feedback  of  the  form  u  =  uxxxxz  4~  uxxyx2y- f  uxyyxy 2  4- 

ui /yy  V3  • 
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5.2.  The  case  bi  =  0 


Next,  we  consider  the  case  in  which  Hypothesis  2B  holds,  i.e.,  by  =  0  and  An  is  as  in 
(38).  Let  the  control  input  be  of  the  form 


u  =knx  +  k12y  +  I<2(  +  U(x,y,£), 


(51) 


where  U  is  defined  in  (45). 

Assume  that  A22  +  b2K2  is  stable.  From  Section  3,  the  stability  of  the  origin  of  (9) 
agrees  with  the  stability  of  the  origin  of  the  reduced  model 

x  =  Slxy  +  f(x,y,Eix  +  E2y  +  h(x,y))  (52) 

V  =  £l2x  +  g(x,y,Eix  +  E2y  +  h(x,y)).  (53) 

Here,  E  =  (EX,E2)  and  h(x,y )  are  the  solutions  of  Eqs.  (15)  and  (17),  respectively,  with 
K\  —  (fen,  k\2). 

Since  (A22-\-b2K2)  is  stable,  matrices  (A22+b2K2)2 +£IiQ2I  and  (A22+b2K2)2 +4:£liQ,2I 
are  then  both  invertible. 

Let 

H(x,y)  :=b2U(x,  y,  E\x  +  E2y)  +  G(x,y,Eyx  +  E2y) 

~  f(x,y,Eyx  +  E2y)E\  -  g(x,y,Exx  +  E2y)E2 
-x2Hxx  +  xyHXy  +  i/Hyy  +  0 ( 1 1 ( x ,  y ) |  | 3 ) .  (54) 

As  before,  we  take  h  to  be  of  the  form  (44).  Solving  Eqs.  (15)  and  (17),  we  have 

Ey  =  —  {(A22  E  b2Ii2)2  +  £li$l2I}  1  {k\i(A22  +  b2I{2)  —  £l2ki2l}b2  (55) 

E2  =  —  {(A22  +  b2Ii2)2  +  1  {^’12(^22  +  b2I(2)  +  0,yknl}b2  (56) 

and 

hxy  — {(A22  +  b2I\2)2  +  4O1O2-O  1  {2Q.2Hyy  —  2QyHxx 


—  (A22  +  b2K2)Hxy}, 

(57) 

=  —  (A22  +  b2K2)  1(HXX  +  £l2hxy): 

(58) 

=  ~  (A22  +  b2K2)~1(Hyy  ~  £llhXy). 

(59) 
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The  reduced  model  (52)-(53)  is  obtained  as 


X  — y  T  fxxx  T  fxyXy  +  fyy9  T  fxx xx 

+  fxxyx  V  T  fxyyxV  T  fyyyV  T  0(\ |(x,  y)\ |  )  (60) 

V  —  02X  T  9xxx  T  §xyxV  T  9yy9  T  9  XXX'E 

T  9xxyx  y  +  gxyyXy  +  9yyy9  T  0{  |  |(.T,  y)  \  |  ).  (61) 


Here,  ,  fijk  and  gijk,  i,j,  k  €  {x,  y,  z },  denote  the  new  versions  of  the  quadratic  terms 

and  cubic  terms,  the  values  of  which  are  given  in  Appendix  A. 

Referring  to  the  stability  criterion  (41)  and  the  preceding  discussions,  we  obtain  sta¬ 
bility  conditions  for  the  control  system  (9)  as  summarized  in  the  following  lemma. 

Lemma  7.  Suppose  Hypothesis  2B  holds  and  that  the  control  input  is  of  the  form  (51) 
with  nonlinear  function  U  as  in  (45).  Then  the  origin  of  Eq.  (9)  is  asymptotically  stable 
if  A-22  +  b2K2  is  stable  and 


gXy(~ 9xx  +  ±§yy)  -  fzyi-^fxx  +  ^ fyy)  +  ^Jxx9xx 

20,2  f  *  j;  02  £  1  02  ^  n 

j^T"  fyydyy  +  3(/xxx  +  0^  Jxyy  T  ~^9xxy  +  Q^9yyy)  <  0. 


(62) 


Remark  2.  From  (62)  and  Appendix  A,  we  observe  that  only  quadratic  terms  of  the 
function  G,  and  the  linear  and  quadratic  terms  of  the  control  input  u  contribute  to  the 
stability  conditions.  A  linear  and/or  quadratic  feedback  stabilizing  control  law  readily 
follows  from  Lemma  7.  Moreover,  a  stability  criterion  for  the  uncontrolled  version  of 
system  (9)  is  also  implied  by  Lemma  7  by  letting  u  =  0. 

Similarly,  although  Lemma  7  addresses  the  design  of  a  linear  feedback  stabilizing 
control  law,  such  a  linear  stabilizing  control  law  need  not  exist.  In  the  next  result,  we 
consider  a  special  case  of  which  the  non-critical  state  £  of  system  (9)  is  a  scalar.  Since  £ 
is  a  scalar,  as  observed  from  Eqs.  (55)-(56),  we  always  have  solutions  for  the  control  gains 
En  and  k12  for  arbitrary  given  values  of  E\ ,  E2  and  K2-  According  to  the  formulations  as 
in  Appendix  A,  we  can  select  E\  =  0  and  E2  large  enough  (or  E2  =  0  and  E\  large  enough) 
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such  that  the  condition  (62)  in  Lemma  7  holds  while  <  0  (/^  <  0).  The  following 
result  hence  readily  implied  by  Lemma  7. 

Corollary  7.  Suppose  the  non-critical  state  £  is  a  scalar  and  Hypothesis  2B  holds.  Then 
there  is  a  purely  linear  feedback  which  asymptotically  stabilizes  the  origin  of  (9)  if  either 

htt  <  0  or  <  °- 

Referring  to  Eqs.  (54)-(56),  for  the  general  case  of  which  the  state  £  of  system  (9) 
may  not  be  a  scalar,  we  have  H(x,y )  =  b2U(x,y,0)  +  G(x,y,  0)  while  kj  i  =  k\2  =  0  and 
I<2  =  0.  A  purely  quadratic  stabilizing  control  law  can  then  be  obtained  as  follows. 

Corollary  8.  Assume  that  Hypothesis  2B  holds,  A22  is  stable  and  the  origin  of  system  (9)  is 
unstable.  Then  a,  purely  qxiadratic  stabilizing  feedback  in  the  form  u  —  uxxx2+uxyxy+uyyy 2 
exists  for  the  origin  of  (9)  if  one  of  the  following  three  conditions  holds: 

(i)  M0A22b2  0,  or 

(ii)  {(3/x£  +  \gyi)A^  +2Q1Mo}b2  ±  0,  or 

(hi)  {2kl2M0  -  +  9ye,)A22}b2  ^  0,  where 

Mo  =-{kl2(2gyz  -  8fxt)A22  +  +  9x{\{A\2  +  4flifl27)  1.  (63) 

We  note  that  Aeyels’  stabilization  conditions  for  a  third-order  system  [1]  are  special 
cases  of  those  given  in  Corollary  8.  Moreover,  Corollary  8  easily  extends  to  the  case  in 
which  A22  is  not  stable  but  the  pair  (A22,  b2)  is  stabilizable.  This  involves  use  of  an 
additional  linear  feedback  Ab£  to  ensure  the  existence  of  a  locally  invariant  manifold  and 
the  stability  of  the  Jacobian  matrix  of  Eq.  (9b).  Furthermore,  we  note  that  the  stabilizing 
control  laws  obtained  in  Corollaries  6  and  8  agree  with  those  obtained  by  Abed  and  Fu  [6], 
where  an  asymptotic  expansion  method  based  on  bifurcation  analysis  is  used  for  controller 
design. 

6.  CONCLUDING  REMARKS 

In  this  paper,  the  center  manifold  reduction  technique  has  been  applied  to  the  smooth 
feedback  stabilization  of  nonlinear  systems  in  two  critical  cases.  The  stabilizing  control 
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laws  were  obtained  in  a  two  step  composite-type  design.  Linear  stability  for  the  noncritical 
state  £  is  first  ensured,  then  the  remaining  control  gains  are  chosen  to  stabilize  the  origin 
of  the  reduced  model  whose  eigenvalues  all  lie  on  the  imaginary  axis.  Stabilizing  control 
laws  have  been  designed  in  linear  and/or  nonlinear  feedback  form. 
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APPENDIX  A 

The  coefficients  in  the  Taylor  expansions  of  f,g  are  given  below  in  terms  of  those  of 
/,  g.  Here,  p  denotes  either  /  or  g,  and  i  j  for  i,j  e  {x,y}  with  £p.]  =  E\,  and  E[y]  =  E2. 

Pa  -  Pa  +  PiiE[i\  +  E{i]P^E[i] 

Pij  -  Pi j  +  Pi£E[j]  +  PjiE[i]  +  2E'[i]PttE[j} 

Pm  —  pm  T  Pn%E[i]  T  T  P^^(.E[i],>  E[i]j  E[i]) 

+  Pi£  ha  +  2E^p^ha 

pnj  —  Pj^hn  ~b  pi^ij  T  ^E[j] p^^ha  T  2E[jj p^^htj 

+  puj  +  Pij$,E[i]  +  pu^E[j]  +  E[qpjtfE[{\ 

+  2E'[i\PittE[j]  +  3PttdE[i\’E[i]’E[j})- 
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